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Abstract. We perform a detailed analytical and numerical study of the dynamics of perturbations (vor¬ 
tex/aperiodic mode, Rossby and spiral-density waves) in 2D compressible disks with a Keplerian law of rotation. 
We draw attention to the process of spiral-density wave generation from vortices, discussing, in particular, the 
initial, most peculiar stages of wave emission. We show that the linear phenomenon of wave generation by vortices 
in smooth (without inflection points) shear flows found by using the so-called non-modal approach, is directly 
applicable to the present case. After an analytical non-modal description of the physics and characteristics of the 
spiral-density wave generation/propagation in the local shearing-sheet model, we follow the process of wave gen¬ 
eration by small amplitude coherent circular vortex structures, by direct global numerical simulation, describing 
the main features of the generated waves. 
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1. Introduction 

The dynamics of vortices in astrophysical disks has re¬ 
cently received much interest both because vortices in 
protoplanetary disks can represent aggregation regions of 
solid particles for the eventual formation of planets (Barge 
& Sommeria 1995) and more generally for understanding 
accretion disk dynamics and the basic problem of angu¬ 
lar momentum transport (Lovelace et al. 1999, Li et al. 
2000). Several works have been devoted to the analysis of 
the possibility of forming and maintaining coherent vor¬ 
tex structures in the strongly sheared flow pertaining to a 
Keplerian disk, both in barotropic conhgurations, where 
the initial potential vorticity perturbation is conserved 
(Bracco et al. 1999, Godon & Livio 1999,2000, Davis et al. 
2000, Davis 2002) and in baroclinic situations where one 
can have vorticity generation (Klahr & Bodenheimer 2003, 
Klahr 2004). In the incompressible case it has been shown 
that coherent vortex structures can indeed form (under 
conservation of potential vorticity) and anticyclonic vor¬ 
tices can survive longer than cyclonic ones (Bracco et al. 
1999) and give rise to the appearance of Rossby waves in 
the system (Davis et al. 2000). The effects of compressibil¬ 
ity have not yet been fully analyzed and require rigorous 
study. Godon & Livio (1999) performed two-dimensional 
time-dependent numerical simulations of vortices in vis¬ 
cous compressible Keplerian disks. Vorticity waves are 
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considered as one of the constituents of (anticyclonic) vor¬ 
tex dynamics, but without specification of the wave prop¬ 
erties and any analysis of their genesis and dynamics (the 
subject of the study was the stability and lifetime of vor¬ 
tices). Davis (2002) performed fully compressible numeri¬ 
cal simulations of the dynamics of a single vortical struc¬ 
ture in a Keplerian disk flow and has reported the gen¬ 
eration of outward-moving compressible waves by the co¬ 
herent vortex (the generation was attributed to nonlinear 
processes) pointing out the potential importance of this 
phenomenon for vortex dynamics. Johansen, Andersen 
& Brandenburg (2004) considered the dynamics of non¬ 
linear vortices by numerical 3D simulations in the local 
shearing sheet approximation, also observing indications 
of wave generation. Klahr & Bodenheimer (2003), per¬ 
forming 2D and 3D hydrodynamical simulations of pro¬ 
toplanetary disks, found that a radial entropy gradient 
can generate Rossby waves which eventually break into 
vortices. Klahr (2004), by a linear stability analysis sim¬ 
ilar to the one presented in this paper, in fact showed 
that a radial entropy gradient leads to continuous gener¬ 
ation of potential vorticity and to a transient swing-like 
amplihcation of the vortical/aperiodic mode, without the 
subsequent decay observed in the linear barotropic config¬ 
uration. 

The aim of this paper is to investigate the linear dy¬ 
namics of initially-imposed vortical perturbations in 2D 
compressible Keplerian disks to understand the phenom- 
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ena of spiral-density and Rossby wave generation and 
to study the initial, most peculiar stages of wave emis¬ 
sion/propagation. In this repect, we want to stress that 
the origins of spiral-density and Rossby waves are abso¬ 
lutely different. 

We base our understanding of spiral-density wave gen¬ 
eration on the papers by Chagelishvili et al. (1997, 2000), 
where a completely linear phenomenon of acoustic wave 
generation by vortex mode perturbation in smooth shear 
flows was for the first time discussed. It was also shown 
that the generated acoustic waves are emitted as two 
symmetric packages in opposite directions from the par¬ 
ent coherent vortex perturbation. This wave generation 
phenomenon is directly applicable to the generation of 
spiral-density waves by vortices in disk flows (Sec. 3). The 
Rossby wave appearance from the vortex mode perturba¬ 
tions, on the contrary, is not connected to any generation 
process. The essence of the Rossby wave formation lies in 
the “breaking of degeneracy” of vortex mode perturba¬ 
tions (Sec. 4). 

Recently, much progress has been made by the hydro- 
dynamic community in the analysis of the linear dynamics 
of smooth shear flows by recognising that the traditional 
normal mode approach does not represent the best tool for 
its study because of the non-normality of the linear oper¬ 
ators in the flows (cf. Reddy, Schmid & Henningson 1993, 
Henningson & Reddy 1994, Gustavsson 1991, Farrell & 
loannou 1993, Craik & Criminale 1986). This means that 
the corresponding eigenfunctions are not orthogonal and 
strongly interfere. The knowledge drawn from the analy¬ 
sis of separate eigenfunctions and eigenvalues is far from 
complete. A correct and full description of the shear flow 
phenomena needs knowledge of the interference processes 
which cannot be easily taken into account in the frame¬ 
work of modal analysis. These circumstances led to the 
development of the so-called non-modal approach that 
could reveal several unexpected phenomena, which were 
overlooked by the normal mode analysis. The non-modal 
analysis is a modification of the initial value problem. It 
involves the change of independent variables from the lab¬ 
oratory to a moving frame and the study of the temporal 
evolution of spatial Fourier harmonics (SFH) of pertur¬ 
bations without any spectral expansion in time. The re¬ 
sulting linear dynamics is much richer than that expected 
from the stability predictions of the modal approach and 
these linear processes may also play a fundamental role 
in the full dynamics of smooth shear flows. For example, 
it has been shown that vortical perturbations can have a 
phase of transient amplification, which is at the base of 
the so-called “bypass” concept of the onset of turbulence 
in planar Couette flows (Baggett et al 1995, Gebhardt & 
Grossman 1994, Farrell & loannou 1993, Henningson & 
Ready 1994, Grossmann 2000, Chagelishvili et al 2002, 
Chapman 2002). Another phenomenon of interest here is 
the generation of waves by a vortex/aperiodic mode of 
perturbations, which has been discussed for the first time, 
as we saw above, by Chagelishvili et al. (1997) and ana¬ 
lyzed by Chagelishvili et al. (2000) and Farrell & loannou 


(2000) in the simplest planar configuration. In this paper 
we extend the non-modal analysis of this process to the 
case of 2D Keplerian disk flow (that is a natural example 
of smooth shear flows) in the shearing sheet approxima¬ 
tion, giving the theoretical basis to understand the results 
of simulations like those presented by Davis (2002). We 
also perform numerical simulations of the vortex dynam¬ 
ics in Keplerian disks: we consider a small amplitude (to 
include only linear processes), coherent circular vortex and 
focus on wave generation and emission, on the properties 
of density and vorticity fields and on the description of the 
wave propagation trajectories. Physical model and equa¬ 
tions are presented in Sec. 2. The non-modal analysis of 
spiral-density wave generation/emission by vortex mode 
perturbation is presented in Sec. 3. The appearance of the 
Rossby wave from vortex mode perturbation is outlined in 
Sec. 4. Numerical simulations of wave emergence from a 
coherent circular vortex perturbation is presented in Sec. 
5 and in the last section we summarize our results. In 
the Appendix we describe the method for the selection of 
initial values of pure vortex mode perturbations for our 
numerical study. 


2. Physical model and equations 


In cylindrical coordinates, the basic equations will read as 
follows: 
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We consider a two dimensional (2D) inviscid, differen¬ 
tially rotating flow. The flow has axial symmetry and a 
rotation axis parallel to the z axis: Vq = (0, Vo,^,0), 
with = rn(r). In this flow configuration the basic 

force balance is described by: 


rVfl{r) = 


1 dPp i94>(r) 
po dr dr ’ 


( 6 ) 


We choose a distribution of angular velocity that follows 
the Keplerian law 


n^(r) (x ^ 


3 ’ 


( 7 ) 


and a gravitational potential that balances the centrifugal 
term 
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the pressure is then constant and equal to Pq. 

Our 2D model of a Keplerian disk retains the dynam¬ 
ical effects of differential rotation, while discarding geo¬ 
metrical and thermodynamic { Pq-, Pq = const) compli¬ 
cations. The sound speed is defined by: 

Cs = = const. (9) 

Po 

with the specific heat ratio 7 = 5/3. Even though our 
flow is two-dimensional and has no z dependence, it is 
useful for the discussion to introduce the parameter 

H = cs/n , (10) 


we derive, from system HE the equations governing the 
linear dynamics of perturbations in the local Cartesian 
frame: 


d „ . 9 
m + 


Vy 2(flQ -t- At)v^ -\- 


Po dy 


= 0, (18) 


d d\ , dv' 

^+2^^^jp+Po + 


(19) 


that represents the equivalent thickness of a thin disk 
model. 

In a Keplerian disk the basic flow vorticity is (cf. 
Tagger 2001): 



APf + 2D 



D 

' 2 ' 


( 11 ) 


Hence, following the standard method of non-modal 
analysis (cf. Goldreich & Lynden-Bell 1965; Goldreich & 
Tremaine 1978), we introduce the spatial Fourier harmon¬ 
ics (SFH) of perturbations with time-dependent phases: 


Thus, the mean vorticity gradient across the disk is 
nonzero and should not a priori be neglected. 

3. Linear dynamics of vortex mode - Generation 
of spiral-density waves 


f Vx'{r,t) ' 


" Vx{kx{t),kyA) ' 

1 Vy'iyA) 

L 

^y{kx(t)^kyA) 

1 p'(r,t) 

" 

p{kx{t),ky,t) 

1 P'(Di) , 


^ g{kx{t),ky,t) ^ 


exp(ifc2,(t)a; -I- ikyp) 


( 21 ) 


Here we study the linear dynamics of 2D small-scale per¬ 
turbations (i.e. with characteristic scales in radial and az¬ 
imuthal directions much less than the radial character¬ 
istic scale of the background/disk flow) in the shearing 
sheet approximation (e. g., Goldreich & Lynden-Bell 1965; 
Goldreich & Tremaine 1978; Nakagawa & Sekiya 1992). In 
this case the dynamical equations are written in the local 
co-moving Gartesian co-ordinate system: 

x = r-ro; y = ro(<^ - Dot); -,-^< 1 , ( 12 ) 

ro ro 

where (r, </) are standard cylindrical co-ordinates and Dq 
is the local rotation angular velocity at r = rp: 

D(r) = Do-I-^(r - ro) = Do-I-2H—, (13) 

or ro 

where 


is the standard Oort constant which describes the mean 
velocity shear parameter in the local frame. In the case of 
a Keplerian rotation law: 

H = -^Do<0. (15) 

The main assumption of the shearing sheet approximation 
is the neglect of the basic flow vorticity gradient; the only 
gradient being considered is differential rotation. 
Introducing linear perturbations 


kx(t) = kx(0) - 2Akyt = kx(0) + ^Dofcyt. (22) 

The streamwise/azimuthal wavenumber remains con¬ 
stant and the streamcross/radial wavenumber changes at 
a constant rate; in the linear approximation SFH “drifts” 
in the k-plane (wavenumber plane). The effect of this 
change is that lines of constant phase (wave crests) are 
rotated by the basic flow. 

Substitution of Eqs. Tmm into Eqs. mm yields the 
system of ordinary differential equations that governs the 
linear dynamics of SFH of perturbations in the described 
flow: 

= 2noVy (t) - ikx (t) —, (23) 

dt Po 

±M=-2{no + A)vx{t)-iky^, (24) 

at " Po 

+ ipo [kx{t)vx{t) + kyVy(t)] = 0, (25) 

P{t) = clgit). (26) 

Introducing D{t) = ig{t)/po, Eqs. (23-26) become: 

= 2DoWy(t) - kx(t)c^D(t), (27) 

= _2(Do + A)vxit) - kycjDit), (28) 

dD(t) 


dn 


(14) 


V = Vo-bv', P = Po+p\ p = po + p\ 


(16) 


dt 


kx{t)Vx{t) -I- kyVy{t). 


(29) 
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Klahr (2004) obtained a similar system of equations for a 
configuration with a radial entropy gradient in the basic 
state, in that case he gets a fourth order system, while 
in our case, as we will see, the system can be reduced to 
second order. 

The latter system is characterized by the important 
time invariant: 

kx{t)vy{t) - kyVxit) + 2(^0 + A)D{t) = X (30) 

that follows (for SFH of perturbations) from the con¬ 
servation of potential vorticity. This time invariant, in 
turn, indicates the existence of the vortex/aperiodic mode 
in the perturbation spectrum whose dynamics represent 
the main interest of this Section. A radial entropy gradi¬ 
ent would introduce a source term for potential vorticity 
(Klahr 2004) and the quantity above would no longer be 
invariant. 

We can define the spectral density of the total energy 
in the k-plane as: 

£k(() = f (d+»^) + ^c^ (31) 


to X and goes to zero when X = 0. The existence of 
the vortex/aperiodic mode, as discussed in more detail in 
Section 4, is the result of the main assumption of the shear¬ 
ing sheet approximation, i.e. the neglect of the basic flow 
vorticity gradient. Also, the vortex mode acquires nonzero 
divergence and nonzero density perturbation ( ^ q, 

where is the density perturbation of the vortex mode 

as defined below). 

Thus, we can decompose any perturbation as the sum 
of the wave and vortex mode perturbations: 

Vx = ; Vy = . D = + JjM . 

(36) 

Clearly, the character of the dynamics depends on 
which mode of perturbation is imposed initially in Eqs. 
mm- pure waves (without mixture of aperiodic vor¬ 
tices), pure aperiodic vortices (without mixture of waves) 
or a mixture of the two. 

Here we present a numerical and qualitative study of 
the linear dynamics of an initially excited SFH of vortex 
mode perturbation. The dynamics is defined by the value 
of the parameter: 


where the two terms correspond to the kinetic and poten¬ 
tial energies of SFH. 

The numerical study of SFH dynamics is governed by 
Eqs. (I27I30II . However, for the fundamental comprehen¬ 
sion of the physical nature of the perturbations and their 
linear dynamics in the flow, we rewrite them in another 
form: from Eqs. (™nii one can get the following second 
order inhomogeneous differential equation for Vy{t): 

(0 = kx {t)clX, (32) 

and the equations expressing Vx{t) and I?(t) as functions 
of Vy(t) and its time derivative: 




4(Oo -l- A)^ 


kyC^S 


kx{t)kyclvy{t) - 2(Ho + A) 


dvyjt) 

dt 


kyCgX 


(33) 


Dit) = A) “ kyVx{t) - X] , (34) 

where 

uygoik) = 4f2o(Ho -I- A) -I- fc^(t)cg = fig -I- , (35) 

and k‘^{t) = kl{t) + ky. 

Eg. 1321 describes two different modes of perturbations: 

(a) Spiral-density wave mode ( r'y'^^), that is described 
by the general solution of the corresponding homogeneous 
equation. Notice that the frequency of the wave uj{t) is 
time-dependent. 

(b) Vortex mode ( that is aperiodic, originates from 

the equation inhomogeneity { kx{t)c^X) and is associated 
with the particular solution of the inhomogeneous equa¬ 
tion. The amplitude of the vortex mode is proportional 


R{t) = 


2|A| 

4Ho(Ho -I- A) -I- [kx{t) -I- ky]cg 

3 

2^1 +[km+ kl]H^ ’ 


(37) 


where H is the equivalent thickness of a thin disk model 
defined in Eq. (trnil . R{t) is a time-dependent parameter 
and consequently, the character of the dynamics is also 
time-dependent. The final result of the dynamics is defined 
by the maximum value of R{t) (achieved at t = t* , when 

kx{e)=0): 


R 


max 


3 

2^1 +km 


3 

2v'l -b {jnH/RY ’ 


(38) 


where m is the azimuthal wave-number of perturbation 
SFH. 

We carried out numerical integration of Eqs. m 
EHl for an initially imposed leading vortex mode SFH 
( kx(0)/ky < 0), without mixture of waves. Initial vortex 
SFH perturbations are defined by the equations presented 
in the Appendix. The calculations show a fundamental dif¬ 
ference in the dynamics of vortices at low and moderate 
Rrriax - 

At Rmax ^ 1 all dynamics occurs at small R{t), 
at low shear rates compressibility is negligible and only 
one phenomenon occurs, i.e. transient growth of SFH of 
aperiodic vortices. The leading SFH gains energy from 
the mean flow and increases its amplitude; when the 
SFH becomes trailing, it gives back the energy to the 
flow and decreases its amplitude. It is obvious that the 
transient growth of vortex perturbations exists in the 
spectrally/asymptotically stable flow system (the rotation 
curve derived from Kepler’s law is extremely stable based 
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Fig. 1. The evolution of SFH normalized velocity and den¬ 
sity perturbations {vx/cg, Vy/cg and D = ig/po) and its 
normalized energy ( ifk(t)/£^k(0) ) at Rmax = 1-06 
(m = R/H in terms of the azimuthal wavenumber). In 
the dynamics of the initially imposed vortex SFH (with 
kx{0)/ky = —15, ky = 1) transient growth is dominant, 
but the wave generation starts to be noticeable. 





Fig. 2. The evolution of SFH normalized velocity and den¬ 
sity perturbations {vx/cg, Vyjcg, D) and its normalized 
energy ( i?k(t)/ifk(0) ) at Rmax = 1-34 {m = 0.5R/H 
in terms of the azimuthal wavenumber). Evolving in the 
shear flow, the initially imposed vortex mode SFH (with 
kxi0)/ky = —15, ky = 0.5 ), gaining energy from the mean 
flow and amplifying, retains its aperiodic nature until 
kx{t)/ky < 0. After time t* = 10 (at which fcy(lO) = 0), 
we observe the appearance of the wave SFH. 


on the Rayleigh criterion). The transient growth of two 
dimensional vortical/aperiodic perturbations in an incom¬ 
pressible Keplerian flow was found by Lominadze et al. 
(1988). Later on, the existence of transient growth was 
confirmed by Fridman (1989) and loannou & Kakouris 
(2001). The importance of transient growth for the on¬ 
set of hydrodynamic turbulence in Keplerian disks is dis¬ 
cussed by Chagelishvili et al. (2003) and Tevzadze et al. 
(2003). 

At moderate/large Rmax, the dynamics are richer: 
our calculations show the appearance of a conversion of 
vortices to spiral-density waves at moderate shear rates. In 
Figs. 1-2 we present the results of calculations at Rmax = 
1.06; 1.34 ( respectively m = R/H and m = 0.5R/H 
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Fig. 3. The dynamics of 


vortex (wi''^ / Cg , Vy'''^ / c. 


yA 


10 20 


D, their 


^x/Cg, 

and wave parts 
( Vy'"'’/cg, for a pure vortex initial per¬ 

turbation ( uA^(O), '(;y'^^(0), Zl("'^(0) = 0) with the same 
parameters as in Fig. 2. The process may be interpreted 
as an abrupt emergence of wave SFH from the related vor¬ 
tex SFH at t = t*. The amplitude of the generated wave 
v^\t*) smoothes the jump in the aperiodic mode (see 
Eq. 39). 


in terms of the azimuthal wavenumber) . We can start to 
notice the conversion at about Rmax—^ (i-e., at about 
kyH ~ 0.7), but the overall behavior is still dominated by 
the vortex transient growth, as is evident from the lower 
right panel in Fig. |3 that shows the SFH energy growth 
and decrease. The process already becomes dominant in 
the SFH dynamics at Rmax > 1-2 (i.e., at about kyH < 
0.5). 

In Fig. 2 we present the evolution of non-dimensional 
perturbed quantities ( Vx/cs,Vy/cs, D ) and of the nor¬ 
malized energy ( ifk/i?k(0) ) of a SFH, for the case with 
Rmax = 1-34 . In this case, at the beginning R{t) ^ 1 (as 
kx{0)/ky = —15, i.e., {kx{Q)/kyY » 1), and the initial 
stage of the evolution is incompressible. Evolving in the 
shear flow, the vortex SFH gains energy from the mean 
flow and amplifies, while retaining its aperiodic nature. 
R(t) increases, becomes moderate and compressibility 
comes into effect. As a result, we observe the appearance of 
an oscillating part of SFH, i.e. we observe the appearance 
of spiral-density waves. Thus, the linear dynamics of a 
vortex mode perturbation is followed by the generation of 
a spiral-density wave. When R{t) is moderate, the time 
scales of the vortex and wave SFH are comparable and 
the perturbations are not separable/distinguishable, i.e. 
we have a mix of SFH of aperiodic and oscillating modes. 
Subsequently, R{t) becomes small again, the time scale 
of the wave SFH becomes much shorter than that of the 
















































6 


G. Bodo et al.: Spiral density wave generation by vortices in Keplerian flows 


vortex SFH and the modes become clearly distinguishable 
again. 

In spite of what was said above, we carry out a sepa¬ 
ration of the modes also in the vicinity of time t*, when 
kx{t*) = 0 and R{t) is not small. The performed separa¬ 
tion allows us to establish the initial characteristics of the 
generated wave SFH, that, in turn, determine the further 
dynamics/energetic of SFH, even quantitatively. 

In Fig. Owe present separately the dynamics of Vyjcs, 
Vxjcs-, D and of their vortex /cs, /cs, 
and wave components /cg, 

for a pure vortex/aperiodic initial perturbation 
( (0) = 0), with the same pa¬ 

rameters used for the case in Fig. 2. From Fig. 0 we 
notics: 

- is an odd function of the argument t — t*; 

- and are even functions of the argument 

t - t*; 

- At time t = t* , at which kx{t*) = 0, the vortex mode 
abruptly gives rise to the corresponding wave SFH (see 
Fig. 3, top panels); 

- At the moment of wave SFH emergence, the phase is 
such that the value of v^\t ^ t*+), D{w){t t*+) = 0 
and v^\t t*+) is maximal (see Figs. O bottom 
panels). More precisely, 

|uW(t ^ t*+)\ = 2|uW(f ^ t*-)\ = 2\vy{t*)\ ; (39) 

- The excited wave harmonic (satisfying condition 
kx{t)/ky > 0 ) gains energy from the shear flow and is 
amplified; 

- the vortex SFH after the transient growth gives energy 
back to the background flow. 

This process can be altered by the action of nonlinear 
forces, or by the presence of a radial entropy gradient 
in the Keplerian flow, in which case vorticity is contin¬ 
uously generated, feeding the vortex with Rossby waves 
(see Klahr 2004). 

4. Qualitative study of Rossby waves 

In the previous section we found that the compressible 
shearing sheet approximation results in a third order sys¬ 
tem of ordinary differential equations (see Eqs. I27I291 1. 
or, equivalently, in a second order inhomogeneous dif¬ 
ferential equation (see Eq. E2Ii- This reduction is con¬ 
nected to the conservation of potential vorticity and re¬ 
sults in the fact that one of the perturbation modes is ape¬ 
riodic/vortex. Consequently, the approximation involves 
two spiral-density waves propagating in opposite direc¬ 
tions and one vortex mode. The vortex mode and spiral- 
density waves are no longer independent: they are coupled 
by the strong differential character of Keplerian rotation. 
As is shown in the previous section, the vortex mode un¬ 
conditionally generates spiral-density waves, but the op¬ 
posite does not occur, i.e. the waves do not a generate 
vortex/aperiodic mode. 


A question arises: how does the picture change when 
the shearing sheet approximation is not valid, especially, 
when the mean vorticity gradient across the disk is at 
work? We consider this problem qualitatively, borrowing 
some estimates from Tagger (2001). At first, the linear 
dynamics of perturbations is described by a third order 
system of ordinary differential equations that, however, 
cannot be reduced to a second order inhomogeneous dif¬ 
ferential equation, as in the previous case. In this case all 
the perturbation modes are waves: two spiral-density and 
one Rossby wave. The frequency of the first ones is given 
by Eq. 1231 We can get the frequency of the Rossby wave 
from Eq. (16) of Tagger (2001): 

^ ^ 3 H ky 

’ kl{t)+kldr ArkUt) + ky ^ ’ 

From Eq. (Bnil it follows that in a thin {H -C r) 
Keplerian disk, for perturbations with characteristic sizes 
of the order of the disk thickness {kx,ky ~ 1/F?), the 
Rossby wave frequency is much less than H. Considering 
a short time evolution (~ l/^)i the disk thickness 

size Rossby waves represent a non-oscillating and non¬ 
propagating vortical perturbation. In practice, we have 
a degeneration of Rossby waves to vortex mode perturba¬ 
tions and for them the mean flow vorticity gradient can 
be neglected. 

We present a simple sketch of the “breaking of the de¬ 
generacy” in the wavenumber plane {kx,ky) (see Fig. 4). 
The linear dynamics of a vortical perturbation may be 
described by following each of its SFH in the wavenum¬ 
ber plane. We single out a SFH that, initially, is located 
at point 1 in Fig. 4, for which kx/ky <0; ky < 1/r; 
\kx\ ^ !/?'■) be. a tightly leading vortical perturbation 
with large azimuthal size. Formally, the perturbation rep¬ 
resents a Rossby wave SFH. However, according to Eq. 
(Uni, its frequency is substantially smaller than H. Thus, 
the SFH is aperiodic, i.e., it is degenerated to a vortex 
mode SFH. According to Eq. 123), as kx{t) varies in time, 
the SFH drifts in the direction marked by the arrows (we 
present the drift of the SFH only in the upper half-plane 
ky > 0; since the perturbation is real, there is a counter¬ 
part in the lower half-plane.) Initially, as \kx{t)\ decreases, 
the energy of the SFH grows (see Figs. 1-3). Then, the SFH 
reaches the circled domain where \ky\/ {k^ + ky) < r (point 

2) and its frequency becomes no longer negligible (the 
boundaries of the domain, where breaking of degeneracy 
occurs, are indistinct, but we fixed the domain to clarify 
the analysis). Thus, passing this domain, the SFH acquires 
a wave nature, i.e., the degeneration to a vortex mode is 
broken. At the same time, the Rossby wave SFH contin¬ 
ues its growth that lasts until it crosses the line kx = 0 
(point 3). It then generates the related SFH of a spiral- 
density wave (as it does the vortex mode SFH - see Sec. 

3) . The drift after passing the line = 0 is shown by the 
double line arrows to stress the appearance of the spiral- 
density wave SFH and its further drift equal to the Rossby 
wave drift. At point 4, the Rossby wave SFH leaves the cir¬ 
cled domain and becomes tightly trailing, again degenerat- 
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Fig. 4. Simple sketch of the relation between the Rossby 
wave and vortex mode perturbations in the wavenumber 
plane {K^, Ky), see details in the text. 


ing into an aperiodic/vortex SFH. Thus, the Rossby wave 
appearance from a vortex mode perturbation is not con¬ 
nected to a generation process. The essence of the Rossby 
wave formation lies just in the “breaking of degeneracy” 
of the vortex mode perturbation. If we now choose a vor¬ 
tex mode SFH with ky > Ijr (e.g., at point 1'), it is 
not difficult to understand that the drifting SFH misses 
the circled domain and the breaking of degeneration never 
occurs, but, in crossing the line = 0 (see point 3'), it 
generates its related spiral density wave SFH. 

A supposition about the coupling of Rossby and spiral- 
density waves is expressed in Tagger (2001): “density 
waves generating Rossby waves as they propagate, and 
reciprocally Rossby waves spawning density waves as 
they are sheared by differential rotation”. Chagelishvili 
& Chkhetiani (1995) investigated the linear dynamics of 
planetary atmospheric Rossby waves in the /3-plane ap¬ 
proximation, in a case of zonal flow. They found and de¬ 
scribed the generation of spiral-density waves by Rossby 
waves, but not the opposite. It may be that this asym¬ 
metry of the wave coupling is connected to the fact that 
in the system, we have only one Rossby and two spiral- 
density waves (propagating in opposite directions). The 
same asymmetry exists in case of wave and vortex mode 
coupling (see Sec. 3 and Chagelishvili et al. 1997). In this 
context it would be particularly interesting to consider 
the case where a radial entropy gradient is present (Klahr 
2004), since in this case there are two independent Rossby 
waves. The problem of Rossby wave generation by spiral- 
density waves in a given disk model (containing both an¬ 


gular velocity and vorticity gradients) will be the subject 
of a separate investigation. 

5. Generation of spiral-density waves by coherent 
circular vortex structure - Numerical 
simulations 

5.1. Numerical setup 

In this Section we present 2D numerical simulations of 
Eqs. (I1I5II in case of an initial coherent circular vortex 
perturbation superimposed on an equilibrium Keplerian 
flow. 

Since the emphasis of our numerical simulations is not 
only on the vortex dynamics itself but also on the process 
of wave generation, we need initial conditions correspond¬ 
ing to the vortex perturbations in a pure form. 

One of the ways to select the initial conditions for the 
vortex would be to employ an equilibrium vortex con¬ 
figuration. A class of widely used exact steady analytic 
solutions for single vortex patches is given by Saffman 
(1992). These solutions were also used for the construc¬ 
tion of an elliptic vortex in Keplerian shear flow in nu¬ 
merical simulations by Chavanis (2000). Another class of 
time-dependent vortex solutions in uniform shear flows is 
given by Kida (1981). Stability analysis of the Kida vortex 
revealed some aspects of its non-modal behavior - tran¬ 
sient growth due to the background velocity shear (see 
Meacham et al. 1990). However, these vortex configura¬ 
tions are solutions for incompressible flows and cannot be 
directly used in the compressible setup initially preserving 
the no-wave condition. 

An exact equilibrium vortex solution in Keplerian 
shear flow is given by Goodman et al. (1987). However, the 
specificity of this solution is that it has a non-divergent ve¬ 
locity field, yielding a particular enthalpy distribution to 
balance the compressibility effects. In this sense, this vor¬ 
tex is not purely kinematic by nature, but involves ther¬ 
modynamic forces in the equilibrium configuration. 

Our aim is to perform numerical simulations of vortex 
dynamics for a direct comparison and verification of re¬ 
sults obtained within the non-modal analysis in the wave- 
number space. For this purpose, we need an initial vortex 
configuration which will persist in a medium with constant 
sound speed - a property which is not valid for the equi¬ 
librium vortex derived by Goodman et al. (1987). Hence, 
we provide the algorithm for the construction of a pure 
vortex structure in compressible shear flow. 

For simplicity we chose the initial vortex to be circular 
in shape having in mind that it will be sheared into an 
elliptical configuration during the evolution in the shear 
flow. The compressibility effects on vortex perturbations 
in flows with moderate shear parameters (as in Keplerian 
disks) are appreciable at \kx/ky\ < 1, while the circular 
vortex we want to construct contains SFH from this do¬ 
main of the wave-number space. Therefore, in this case 
we need a refined procedure to construct pressure, density 
and velocity fields that correspond to the vortex mode per- 
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turbations in a pure form. For this purpose we use “pre- 
initial” conditions for the perturbations in as follows: 

( 2 I 2 \ ^ 

- —exp[-(a:^-bj/^)/a^] , (41) 

J 

/ 2 _i_ 2 \ ^ 

v'y{x,y,Q) = \ a?^ ) > (^2) 

p'(cr,2/,0) =0 , (43) 

where defines the amplitude of the perturbation, a its 
size, n the shape of the vorticity distribution and the up¬ 
per and lower signs correspond respectively to anticyclonic 
and cyclonic rotation. Although the kinematics described 
by the above equations seems a pure vortex, this is not 
the case because the absence of the density disturbance 
denotes that the equations above describe a mixture of 
spiral density wave and vortical perturbations, with den¬ 
sity disturbances fully compensating each other. We use 
them to describe a “pre-initial” condition, in the sense 
that through them we defined the distribution of I{kx, ky) 
in wavenumber plane. Then, on the basis of the procedure 
described in the Appendix, we constructed our actual ini¬ 
tial condition, that is an almost pure (with minimal ad¬ 
mixes of spiral-density waves) circular vortex with its cor¬ 
responding density field. The velocity field of our initial 
disturbance has somewhat changed with respect to the 
“pre-initial” one. Our initial vortical perturbation repre¬ 
sents a rich ensemble of leading {kx/ky < 0) and trailing 
{kx/ky > 0) harmonics, but only the presence of leading 
SFH sets the conditions for the generation of spiral-density 
waves. 

We have then superimposed this perturbation on the 
equilibrium Keplerian flow and have performed the simu¬ 
lations using our implementation of the PPM scheme in 
its Eulerian version (Woodward & Colella 1984, Mignone 
et al. 2004). The main characteristics of the scheme are a 
parabolic reconstruction that gives third order spatial ac¬ 
curacy in smooth regions of the flow, a nonlinear Riemann 
solver and time advance based on a characteristic projec¬ 
tion that is second order accurate. We have not used an 
angular momentum conserving form of the (f> component 
of the momentum equation, however, we have tested that 
the total angular momentum is conserved by our code dur¬ 
ing the simulations, with a relative accuracy of 10“® . For 
long-term simulations, e.g. hundreds of orbits, one would 
have to use a scheme that conserves the local angular mo¬ 
mentum density (Kley 1998), but, as we only study the 
evolution of 2.5 orbits and only impose a linear pertur¬ 
bation, it is possible to neglect the precise angular mo¬ 
mentum conservation. Since in this problem there are no 
natural units of length, we have arbitrarily chosen our unit 
of length as the radius at which the vortex center is lo¬ 
cated. The unit of velocity is then chosen as the Keplerian 
velocity Vkv computed at this radius, i.e. at r = 1, and the 
disk density is taken as units of density. In these units the 
computational domain covers the region 0.056 < r < 1.94 


and 0 < 0 < 27r, in polar coordinates, with a uniform 
(both in r and in 9) grid of 1024 x 1024 points. We had 
to use such a number of grid points in order to have a 
good resolution on the vortex also covering a large area 
of the disk. In this way, for example, in the case with 
a = 0.1, we cover the vortex with about 100 points. The 
center of the vortex perturbation is then located at r = 1 
and 6 = 7r/4. In our units, the rotation period at the 
location of the vortex is therefore T = 2tt and our sim¬ 
ulations reach time t = 2.5T. Boundary conditions are 
periodic in the azimuthal direction and outflow in the ra¬ 
dial direction. More precisely, in the radial direction we 
impose the zero-derivative of all the quantities, except the 
azimuthal velocity, for which we continue the Keplerian 
profile and impose the zero-derivative on the perturba¬ 
tion. These conditions do not ensure perfect transmission 
of perturbations impinging on the boundaries, however the 
portions of perturbation that are reflected at the bound¬ 
aries become leading waves that are deamplified by the 
shear and can influence only a small region close to the 
boundaries. The main parameters of our system are then 
Cs, a, €r and n (see Eqs. 41-43 above). We have considered 
two values of the sound speed, two sizes of the vortex and 
two values for n. In Table ^ we report the values of the 
parameters used in our simulations together with some 
other nondimensional parameters that may help in com¬ 
paring the importance of different effects on the vortex 
dynamics. For comparison, in the table we reported the 
same parameters for the simulation performed by Godon 
& Livio (1999) and Davis (2002). 

The second column in Table Ogives the value of the 
sound speed (in units of the Keplerian velocity Vk^ at the 
position of the center of the vortex) and, according to Eq. 
m this is equivalent to the values of the ratio H/r be¬ 
tween the disk thickness and the radius. Case B clearly 
represents the case with the lowest value of H/r, however 
all the cases help to gain more insight into the dependence 
of the vortex dynamics on the parameters: a comparison 
between Case A and Case B allows us to investigate the ef¬ 
fect of the sound speed, case C allows us to investigate the 
effect of the vortex size and case D allows us to investigate 
the effect of a different vorticity distribution ( a value of 
n different from 0 gives a vorticity ring, while n = 0 gives 
a filled distribution). Unlike what we said above, in case 
C we made use of 1096 x 4192 grid points covering the 
domain 0.5 < r < 1.5 and 0 < 0 < 7r/2 and the final inte¬ 
gration time is t = 3.5. With this choice we could keep the 
same resolution on the vortex that we have for the other 
cases. The amplitude of the perturbation in all cases, ex¬ 
cept case C, is e = 0.004, and in case C it is ten times 
smaller, in order to minimize the action of nonlinear terms 
of Eqs. (1-4), since in this paper we want to address the 
linear dynamics of the perturbation. We do this to prove 
the linear character of the spiral-density wave generation 
process (the main aim of our study) and to visualize the 
wave generation and propagation processes clearly. In the 
linear case there is no difference between the dynamics 
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Table 1. Parameters of the simulations and estimated values of other nondimensional parameters at the initial vortex 
position. For comparison we also give values of some parameters for Godon & Livio (1999) and Davis (2002) simulations 



Cs/vkv (= H/r) 

a 

n 

Ad/r 

Ad/H 

Vm/Vkbp 

N 

A 

0.26 

0.1 

0 

0.25 

1 

0.0034 

0.054 

B 

0.08 

0.1 

0 

0.25 

3.1 

0.0034 

0.054 

C 

0.26 

0.01 

0 

0.025 

0.1 

0.00034 

0.054 

D 

0.26 

0.1 

2 

0.5 

2 

0.0015 

0.015 

Godon & Livio 1999 

0.15 

- 

- 

0.05 

0.33 

0.03 

2.4 

Davis 2002 

0.075 

- 

- 

0.25 

3.3 

0.08 

1.3 


of cyclonic and anticyclonic vortices, so, without loss of 
generality we took the vortex rotation to be anticyclonic. 

The parameters Ad/r and Ad/H, as discussed in Sec. 
4, giving the vortex core size Ad in units respectively of 
the radius and of the disk thickness, measure the impor¬ 
tance of the Rossby waves. On the basis of Fn. 1401 one can 
estimate the characteristic oscillation times of the Rossby 
waves in the simulated cases. These times significantly ex¬ 
ceed the time of our simulation - a few revolutions of the 
disk at the vortex location. Consequently, Rossby waves 
do not manifest themselves in our simulations, i.e., the 
real participants of the simulated dynamics are the vortex 
mode and the (generated) spiral-density wave. 

The parameter N in the last column of Table ^ mea¬ 
sures the competition between the distortion by the shear 
flow (linear phenomenon) and nonlinear phenomenona 
that might oppose this distorsion. This competition may 
be estimated by the ratio between the vortex maximum 
velocity vm and the local Keplerian velocity difference 
across the vortex core radius {\dVKep/dr\Ad/2) 

N= __= 4^—. (44) 

\dVKep/dr\Ad/2 Vkbp Ad 

At fV <C 1 the vortex dynamics is completely linear; 
when N becomes of the order of unity or larger, nonlinear 
effects are important. As one can see from the Table, the 
value of this parameter for Godon & Livio (1999) and 
Davis (2002) indicates that in their simulations nonlinear 
effects are important, while all our cases describe strongly 
linear dynamics of vortices. 

5.2. Results 

In Fig. 5 we represent the initial conditions of our simula¬ 
tions (for n = 0), more precisely we show the vorticity (up¬ 
per panel) and density (lower panel) perturbation fields, 
with enlargements of the vortex region. The vorticity field 
of the initial vortical perturbation has axial symmetry, as 
expected, while the density field has a different structure 
resulting from the balance of pressure and Goriolis forces. 
In Figs. 6 and 7, we show the initial phases of the vor¬ 
tex evolution for case A, up to t = 27r, that corresponds 
to one rotation period at the initial vortex location. The 
Figures show respectively images of vorticity and density 
distributions at six different times. We see that the vortex 
is stretched by the background Keplerian flow, but we can 


also observe that, in agreement with the theory presented 
in Sec. 3, two spiral-density waves are being generated 
by the vortex. The stretching is well visualized in Fig.6, 
that shows the vorticity dynamics. The wave generation 
and propagation phenomena are well traced by the two 
arms that are visible on the images of Fig. 7. According 
to what we have shown in previous sections, a vortex SFH 
generates a wave when the SFH crosses the ky axis, i.e. 
when the wavenumber in the local radial direction be¬ 
comes zero (kx = 0, see Fig. 3 and points 3 and 3' on 
Fig. 4). Our vortex is composed of many SFH, and each 
of them, crossing the ky axis, generates the corresponding 
spiral-density wave SFH component. Therefore, since at 
the moment of generation the wavenumber is along the az¬ 
imuthal direction, we can observe the two generated waves 
to be propagating in opposite azimuthal directions. The 
waves, as they propagate, have their wavevectors turned 
by the shear flow towards the radial direction, and in time 
they are carried by the flow and become tightly trailing. 
The appearance of two arms is related to the conserva¬ 
tion of the perturbation action, that initially is zero, since 
the vortex perturbation does not propagate. The total ac¬ 
tion of the generated wave arms should then be zero and 
therefore we have two arms with opposite propagation di¬ 
rections. 

We now define and compute the trajectories of the ex¬ 
cited spiral-density waves by using the group velocity con¬ 
cept (of course, trajectory, as wave group velocity, gains 
in importance only for wave packages when the concept 
of ray is applicable). The trajectories that we will com¬ 
pute will be global, but they will be based on the local 
dispersion relation that can be obtained at each position 
in the shearing-sheet approximation. The validity of the 
trajectories that we will compute has therefore all the lim¬ 
itations of this approximation. Due to this “half-globality” 
we will indicate by kr(t) and kg (instead of kx(t) and ky 
that are used in the local reference frame) the components 
of the wavenumber vector respectively in the local radial 
and azimuthal directions. We rewrite the dispersion rela¬ 
tion of spiral-density waves from Eq. (34) 

LOg(kr{t), kg) = + {k‘^(t) + kg)c^]'^ = 

±[l + H\k^^{t) + kl)]hs (45) 
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Fig. 5. Images of the initial density (upper plot) and vorticity (lower plot) distributions with enlargements of the 
vortex area. 


and then we obtain from it the radial and azimuthal com¬ 
ponents of the group velocity 

_ f 9ujs duJs '\ _ ( _ kr{t)cs _ 

* " [dkrit) ’ dke) ~ [ + HWrit) + k'i] ’ 

i_ kecs 

The two signs identify the two oppositely propagating 
spiral-density waves. We notice that the group velocity is 
time-dependent as kr(t) depends on time. We can now 
compute the trajectory by defining the total wave velocity 
as 

^=Vf=Vf + Vfc (47) 

with the wavenumber that changes according to the equa¬ 
tion 

^ = -2A{r)k. . (48) 

where A{r) is local value of the Oort’s constant. As dis¬ 
cussed in Sect. 3, the group velocity of spiral-density waves 
is directed along the azimuthal direction at the moment of 
their generation, i.e., initially we observe two oppositely 
propagating rays in this direction. We can then follow the 



propagation by numerically integrating Eqs. 14714X11 . In 
Figs. 8 and 9 we plot the trajectories that we have ob¬ 
tained in this way, superimposed on the results of the full 
numerical simulations at time t = 2 for the two cases 
A and C that differ for the initial vortex size. The dif¬ 
ferent trajectories correspond to different values of kg at 
the emission time. The trajectories have a characteristic 
S-shape since the two rays oppositely propagating in the 
azimuthal direction are progressively turned towards the 
radial direction and in the meantime they are carried by 
the background Keplerian flow. We can see that our com¬ 
puted trajectories agree quite well with the results of nu¬ 
merically simulated propagation of spiral-density waves. 
We see that the agreement is better for the case with an 
initially smaller vortex size. In both cases we have part of 
the perturbation that is ahead of the computed trajectory, 
but this part is larger for case A. This partial disagreement 
can be considered as the result of the breakdown of the 
approximations used in computing the trajectories and of 
course, the approximations are better for a smaller size 
vortex. 


In Fig. 10 we again plot trajectories and numerical 
simulation results for case A at a later time {t = 3.75), 
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Fig. 6. Evolution of the vorticity distribution for case A. The six panels show the vorticity distributions at six different 
times. The times are respectively t = 0.31,1,2,3.12,4.5,6.31. At late times the vorticity distribution in the inner 
region becomes quite noisy, however one has to take into account that the calculation of vorticity implies a numerical 
differentiation operation that amplifies the error, larger here due to the grid deformation 


to show that the computed trajectories follow with a very 
good agreement the further perturbation evolution. 

Fig. 7 show a very regular structure of the emitted 
spiral-density waves. Prima facie, they look like shocks, 
but this is not the case. They have a linear origin and 
we can understand the observed regular structure as de¬ 
scribed below. A wide spectrum of wave SFH is generated 
at each moment of time, since the intensity of generation 


is different from zero over a wide range of ky (see Sec. 3), 
and the panels in Fig. 7 represent an integral picture of 
the generated and propagating wave interference. On the 
other hand, according to the graphs of v^\ 
in Fig. 3, the wave harmonics, at the moment of gener¬ 
ation, have very regular phases that are similar to each 
other (specifically, all density perturbation SFH are ex¬ 
cited with zero phase - see graph of in Fig. 3). This 
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Fig. 7. Evolution of the density distribution for case A. The six panels show the density distributions at six different 
times, t = 0.31,1, 2, 3.12,4.5, 6.31 


phase regularity and similarity are the roots of the regular 
interference picture shown in Fig. 7. 

In vortical flows, potential perturbations acquire a vor¬ 
tical nature and this is the case for the spiral-density waves 
that we are considering. We can estimate the wave vortic- 
ity by Eq. (29), rewritten in a “half-global” fashion (as 
was done defining trajectories). For spiral density waves 
we have 2=0: 

kr{t)v^;'\t)-k9vi'^\t) + 2[n{r) + A{r)]D^^'>{t) = 0. (49) 


Defining the wave SFH vorticity as 

= kr{t)vlr^ (t) - kevi'^^ it), (50) 

we finally get: 

=-^nir)D<^'^\t). (51) 

From this equation we can see that the vorticity of the 
propagating wave is determined by the local vorticity of 
the background flow and by the amplitude of the density 
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Fig. 8. Ray trajectories superimposed on the density dis¬ 
tribution for case A at time t = 2. The different rays 
correspond to different values of kg 


perturbation. Both factors grow for the wave propagating 
towards the inner part of the disk and the result is an in¬ 
crease of the wave vorticity in this region. For the outward 
propagating wave, instead, the density amplitude still in¬ 
creases but at a much lower rate than the decrease of the 
background flow vorticity. Therefore the net result will be 
a decrease of its vorticity. This asymmetry is well traced 
in Fig. 6, where we see a vanishing of the outward arms, 
while the inward ones become well pronounced. 

In Fig. 11 we show grey-scale images of the density 
distribution at t = 2 for the cases A,B,C,D. Cases A and 
D relate to equal sizes of parent vortex and sound speeds 
(see Table 1, they differ only for the values of n) and the 
wave propagation pictures are very similar. In case B, due 
to the small value of the sound speed, the wave arms are 
radially compressed and not strongly pronounced. In case 
C, due to the small size of the parent vortex, the generated 
wave SFH has large wavenumbers, their group and phase 
velocities are more or less the same (close to the sound 
speed) and the arms are more uniform (not diffuse). 

Due to linearity, the described phenomena have a short 
life time (a few disk revolutions): in fact vortex distortion 
by the background shear flow cannot be prevented by ve¬ 
locity self-induction (a nonlinear phenomenon). The par¬ 
ent vortex is stretched, becomes tightly trailing (see Fig. 
6) and, in accordance with Sec. 3, is not able to gener¬ 
ate any more spiral-density waves. Also, the trailing vor¬ 
tex mode perturbations {kx{t)/ky > 0) give back energy 
to the flow and disappear (also in accordance with Sec. 
3). The disappearance of the parent vortex is evident in 
Fig. 12, where we show grey-scale images of density at 
later evolution times; the three panels are respectively at 
t = 8.12,11.25,14.5, i.e. up to two and a half revolutions. 


Fig. 9. Ray trajectories superimposed on the density dis¬ 
tribution for case C at time t = 2. The different rays 
correspond to different values of kg. In the image we show 
also a vortex in a symmetric position for better compari¬ 
son. Note that in case C the computational domain covers 
the radial interval 0.5 < r < 1.5. 



-r 0 1 


X 

Fig. 10. Ray trajectories superimposed on the density dis¬ 
tribution for case A at time t = 3.75. The different rays 
correspond to different values of kg. 


Spiral-density tightly trailing waves have extended to the 
whole disk region inside the vortex position, but their gen¬ 
eration seems to have stopped and they have disappeared 
from the outer region. 
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Preliminary results show that nonlinearity makes the 
parent vortex long-lived and the wave generation process 
more permanent (Bodo et al. 2005). Indeed, prevention 
of parent vortex stretching should be the result of posi¬ 
tive nonlinear feedback, i.e. nonlinear regeneration of the 
leading SFH in the vortex spectrum. These replenished 
leading SFH extract the disk flow energy and generate 
continuously spiral-density wave SFH. Consequently, they 
provide a permanent generation of waves. (The concept of 
replenishing of SFH - a root of the positive nonlinear feed¬ 
back - is a theoretical concept, numerical investigation of 
which is under way). 

The results for a single Fourier component, given in 
Section 3, show that the energy of the generated spiral 
density wave increases with time, i.e. trailing spiral den¬ 
sity waves are capable of extracting energy from the mean 
flow. This wave amplification should also be visible in our 
numerical calculations; for this reason, in Fig. 13, we plot 
the perturbation energy (as defined in Eq. 31) vs time. 

In the figure we can see the perturbation normalized 
energy (normalized to the initial value) vs. time for the 
cases A, B and D(case A - solid curve, case B - dot¬ 
ted curve, case D - dashed-dotted curve). The plots go 
up to one revolution time {tmax = 27r), when the action 
of the parent vortex has not yet completely disappeared. 
In this process of wave generation, the vortex acts only 
as a mediator and the waves draw their energy from the 
background flow, both at the beginning, when they are 
generated, and in their propagation, when they are ampli¬ 
fied by the background shear flow. kr{t) of the generated 
wave SFH increases with time and the dynamics of SFH 
becomes adiabatic, therefore the energy of each SFH in¬ 
creases linearly with its frequency: Ak ^ ujsd- Generally, 
according to Eq. 

(52) 

but, at moderate sound speeds (cases A and D), in Eq. (1^^ 
the second term becomes quite rapidly dominant and, as 
time passes, one can write E\^ ~ kr{t)cs, i.e. the energy of 
each SFH increases linearly with time. The increase of k^ 
is determined by the local shear rate, that increases in the 
inner region, at smaller values of the radius r. It is eas¬ 
ily understood that at low sound speeds (case B) in Eq. 

the first term is dominant for a long time and more¬ 
over the increase in kr is slower, since the waves propagate 
more slowly towards the interior and they therefore expe¬ 
rience a smaller shear rate. Consequently, cosd and Ek are 
practically constant. 

A feature of the generated spiral density waves that 
we have discussed above is their regularity, i.e. perturba¬ 
tions associated with the waves have a prominent sign and 
do not average to zero. One consequence of this is that, 
in principle, it may be possible to have a mass flux as¬ 
sociated with them, that would be therefore proportional 
to the wave amplitude. We have tested this possibility by 
plotting the function: 

= f <i0pvr{r,6,t') (53) 
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Fig. 12. Longer time evolution for case A. The three pan¬ 
els show images of the density distributions at three differ¬ 
ent times; from top to bottom the times are respectively 
t = 8.12,11.25,14.5 



as a function of time, at different radial positions (r = 
0.875 and r = 1.125 in Fig. 14, r = 0.1 in Fig. 15). 
This function Q{r,t) represents the total mass that has 
flown through radius r at time t, in units of the mass 
contained inside the initial vortex radial position. Fig. 14 
represents this quantity for positions that are immediately 
inside (r = 0.875 , dashed line) and immediately outside 
(r = 1.125, solid line) the vortex. The figure shows that 
the action of the waves seems to be in the direction of mass 
collecting towards the vortex area. Fig. 15 represents, in- 
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Fig. 11. Comparison of the four different cases at the same time. The four panels show images of the density distribution 
for cases A (top left), B (top right), C (bottom right) and D (bottom left) at t = 2. 


stead, the mass flow induced by the waves in the inner disk 
region (at r = 0.1). Here we see a definite tendency of the 
wave to induce a non-negligible inflow of mass. These plots 
suggest a general trend of the phenomenon, and more re¬ 
alistic results with simulation of longer time dynamics of 
large amplitude coherent vortices are needed (Bodo et al. 
2005). As discussed above, since we are essentially consid¬ 
ering a linear process and since the quantity Al is a first 
order quantity, the values that we get are proportional to 
the initial amplitude of the perturbation we give. The case 
that we present is an anticyclonic vortex; simulations of 
a cyclonic vortex show that all the fluxes are reversed, as 
they should be. Of course, in addition to this mass flux 
that is proportional to the perturbation amplitude, the 
trailing spiral density waves induces an outward angular 
momentum transfer that is quadratic in the perturbation 


amplitude and, therefore, negligible in these linear calcu¬ 
lations. 

6. Discussion and conclusions 

In this paper we present a linear, non-resonant phe¬ 
nomenon of spiral-density wave generation by vortices in 
a Keplerian disk flow that is closely related to the non¬ 
normality of linear dynamics of perturbations in smooth 
(without inflection points) shear flows. This phenomenon 
can be well interpreted by the use of the non-modal ap¬ 
proach, i.e. by following in time the linear dynamics of 
spatial Fourier harmonics of the vortex mode perturba¬ 
tions. We have done this in the shearing sheet model 
and we found that spiral-density wave are generated only 
by the leading vortex mode SFH that meet the condi¬ 
tion kx{0)/ky < 0. The generation takes place at rela- 
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Fig. 13. Energy normalized to the initial value versus time 
for the cases A, B and D (the solid line refers to case A, 
the dotted curve to case B and the dashed curve to case 
D). 



t 

Fig. 14. Plot of the the total mass that has flown through 
radius r at time t versus time, for case A. This figure refers 
to radial positions that are immediately inside (r = 0.875 
, solid line) and immediately outside (r = 1.125, dashed 
line) the vortex. 


tively large wavelengths of SFH - it becomes noticeable 
at about kyH ~ 0.7 and it is dominant in the SFH dy¬ 
namics at kyH < 0.5. At large wavenumbers (small wave¬ 
length) all dynamics of SFH runs at low shear rates and 
only one phenomenon occurs: transient growth of SFH 
of aperiodic vortices. At small wavenumbers (large wave¬ 
length) when kyH < 0.5, the dynamics in addition to the 
transient growth involves conversion of vortices to spiral- 
density waves. After their appearance, the energy of the 
wave SFH increases and the wave SFH itself in time be¬ 
comes tightly trailing kx{t)/ky^l. 

Small amplitude (linear) circular vortex structures 
(that represent a rich ensemble of leading and trailing 


Fig. 15. otal mass that has flown through radius r at time 
t versus time, for case A. This figure refers to r = 0.1. 


SFH) generate two oppositely propagating wave arms in 
accordance with the perturbation action conservation law. 
They have very regular structure and prima facie look like 
shocks, but this is not the case. They have a linear ori¬ 
gin and are the result of the interference of the generated 
spiral-density waves. Their outward and inward propagat¬ 
ing arms are asymmetric due to the disk curved geome¬ 
try and flow asymmetry (H(r) increases outward and de¬ 
creases inward). 

Since spiral-density wave generation is governed by lin¬ 
ear forces, it does not depend on the sign of the vorticity 
of initial aperiodic perturbation and should be equal for 
cyclonic and anticyclonic vortices even in the nonlinear 
regime. Davis (2002) (whose simulations are in a nonlin¬ 
ear regime - see Table 2), reported that a coherent cyclonic 
vortex simultaneously emits a wave as does an anticyclonic 
one (nonlinear cyclonic and anticyclonic vortices differ by 
their lifetime - cyclonic vortices are destroyed in appre¬ 
ciably shorter times than anticyclonic ones). This result 
may be considered as circumstantial evidence for the lin¬ 
ear origin of the wave. 

The fact that the waves are generated by a linear mech¬ 
anism (and not nonlinear ones, as was thought ) apprecia¬ 
bly increases the significance of waves in the disk matter 
inflow/outflow processes and we showed that waves gener¬ 
ated by anticyclonic vortices tend to collect matter in the 
vortex region, which may be important for protoplanetary 
disks, and tend to induce inflow of matter in the inner disk 
regions. 

The described dynamics are short-lived due to the 
small amplitude of the parent vortex. One can specu¬ 
late that the parent vortex could be long-lived and the 
dynamical picture more permanent at larger amplitudes, 
when nonlinear self-induction of the vortex velocity pre¬ 
vents its shear distortion. Indeed, nonlinear regeneration 
of the leading SFH {{kx/ky < 0)) in the vortex spectrum 
can prevent the vortex from stretching. These replenished 
leading SFH gain the disk flow energy and continuously 
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generate spiral-density wave SFH. Thus, permanent non¬ 
linear replenishment of leading SFH of the parent vortex 
could generate permanent disk flow energy extraction and 
wave generation. The effective efficiency of the feedback 
mechanism in permanently maintaing the vortex and over¬ 
coming the linear and nonlinear damping effects inherent 
to a Keplerian disk (Hawley, Balbus and Winters 1999) 
however has still to be proven by numerical simulations 
and analytical means. 

The linear phenomenon of spiral-density wave genera¬ 
tion by vortices investigated in this paper is quite univer¬ 
sal; it occurs in disk flows with high shear rates and should 
be seen in different kind of astrophysical disks (protoplan¬ 
etary disks; quasars; thin and thick galactic and binary 
system disks). 
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7. Appendix - Selection of Initial Values of Vortex 
Mode Perturbations 

The initial values of physical quantities for our two di¬ 
mensional direct numerical simulations of Eqs. (1-4) may 
be composed by the superposition of the background flow 
and local vortex mode perturbation centered at some tq: 

= Vo(r, ^,0) +v'{x,y,0), (A.la) 

P{r,(j>,0) = Po+p'{x,y,0), (A.lb) 

p{r,(l),0) = po + p'{x,y,0). {A.lc) 

where x and y is defined by 

x = r-ro; y = ro{(j) - not); {A.2) 

ro ro 

The background is the Keplerian flow (Vq = 
(0, Vb 0 ,O), = GM/r) with homogeneous pressure 

and density (Pq , Po = const). The rigorous selection of ini¬ 
tial vortex mode perturbation in the (compressible) flow 
requires a refined procedure, as vortex (aperiodic) modes 
acquire nonzero divergence - nonzero density perturbation 
(see discussion in Sec. 3). In order to select initial vortex 
mode perturbations we use the initial value problem in 
the wavenumber plane (k-plane): 


The principal concern with such a solution is its con¬ 
vergence. For instance, this solution diverges at t = t*, 
when kx{t*) = 0 { A ^ 0 ). The area of the diver¬ 

gence strongly depends on the value of the shear rate. 
In our analysis we consider moderate shear rates, when 
A/{kyCs) < I. In this case the above series is conver¬ 
gent for times when \kx{t)/ky\ > 1 . Moreover, the 
first two term of this series appears to be an excellent 
approximation to the exact numerical solution for times 
when \kxit)/ky\ > 2 . Therefore, we can use this ana¬ 
lytic solution to compose the vortex mode perturbations 
that are localized in the wavenumber plane outside the 
\kx{t)/ky\ < 1 area. We calculate the Vy{kx,ky,0) and 
it’s derivative for the vortex mode perturbations in ex¬ 
plicit form: 


'Cyikx^ ky^ 0 ) — 

Cyikx,ky,0) = 


kx SA^kxky 3fc|j -I- 3ky — k^ 

IC^ ^(2 


2Ah 


_ ^2 _ u2 


16 A^kl 


X (^kx 5 ) 5 

(All) 


3kjj + 6k‘jjk‘^, — 2Ak‘ifk‘i + Skf, + 5kt — 2Aklk‘f, 






x'^y 


/C5 


/ <(x,y,0) \ 

v'y{x,y,G) 

p'{x,y,A) 

V p'{x,y,A) ) 


^ Vx)kx , ky , 0) \ 

Vy{kx,ky,0) ^ 

pikx , ky , 0) 

\p{kx,ky,0) J 
X ex.p(ikxX + ikyy). (A.3) 



where 


Cy { kx , ky , 0) 


xX{kx,ky), (A.12) 


^^Vy{kx (t), 



t =0 


Rewriting Eqs. (29,31): 

(t) 

+[^no{no + A) + k'^{t)cl]vy{t) = kx{t)cll, (A.4) 

kx{t)vy{t) — kyVx{t) + 2{no + A)D{t) = const = 2. (A.5) 

Hence, I, Vy{0) and [dvy{t)/dt]^^i^ form the full set of 
initial conditions for Eq. (A.4). We seek the vortex mode 
solution in the following analytic form: 


r;y(0) = no(0)-I- ni(0) -I- ... -|-Wri(0)-|- ... , (A.6) 


where the zero order term is deduced from the stationary 
form of the solution (see Eq. A.4 at A = 0 ) and sub¬ 
sequent terms are derived using the iterative method (see 
also Chagelishvili et al. 1997): 


wo(0) = [fo(0]t=o = 
V„(0) = [Vn{t)]t=0 = - 


kx{t) 


1C 


K. = k'jj + ky + kx, 


Kfj — 


X {^kx ; ky^ , 

(A.7) 

= 0 



Vn-l{t) ' 

dt'^ 

5 

(A. 8 ) 

kl 


(A.9) 

A, 


(A.IO) 


Remaining physical quantities may be deflned using the 
following equations: 


Cxi^X,ky,0) — 


4(Oo + A)2 


y s 


kxkyCgVy{kx , ky , 0) 


2(I7o-f- A)Vy{kx, ky,0) kyCgX{kx, ky,) , 


(A.13) 


p{kx,ky,0) 


Q{kx, ky,0) — 


kxVy{kx, ky, 0)-kyVx{kx, ky,0)-I{kx, ky) . 

(A.14) 

Hence, using any distribution of the potential vortic- 
ity T{kx,ky) (or Uj,(fca;, fcy, 0)) and Eqs. (A.9-14) we can 
construct corresponding initial vortex mode perturbation 
in the wavenumber plane. 


ipo 

2(f2o -|- A) 


where 



